Orthorhombic ABC semiconductors as antiferroelectrics 
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We use a first-principles rational-design approach to identify a previously-unrecognized class of 
antiferroelectric materials in the Pnma MgSrSi structure type. The MgSrSi structure type can be 
described in terms of antipolar distortions of the nonpolar P^/mmc ZrBeSi structure type, and we 
find many members of this structure type are close in energy to the related polar PQsmc LiGaGe 
structure type, which includes many members we predict to be ferroelectric. We highlight known 
ABC combinations in which this energy difference is comparable to the antiferroelectric-ferroelectric 
switching barrier of PbZrOs . We calculate structural parameters and relative energies for all three 
structure types, both for reported and as- yet hypothetical representatives of this class. Our results 
provide guidance for the experimental realization and further investigation of high-performance 
materials suitable for practical applications. 



There has been great progress in the use of first- 
principles methods in the design and discovery of new 
functional materials, most notably for the rapid evalua- 
tion and screening of a large number of known and as- 
yet hypothetical compounds in a target family already 
known to contain compounds exhibiting the desired func- 
tional behavior [US]. A greater challenge is to develop 
ways to identify new functional materials in families in 
which the desired behavior has been previously unrec- 
ognized. In recent work, we used a combined crystallo- 
graphic database/first-principles approach [5] to identify 
semiconducting members of the family of compounds in 
the LiGaGe structure type as a previously unrecognized 
class of ferroelectrics, characterized by spontaneous po- 
larizations and barriers to polarization switching compa- 
rable to the much-studied ferroelectric ABO3 perovskite 
oxides [SJE]. Here, we apply this approach to identify a 
previously unrecognized class of antiferroelectrics, open- 
ing the way to increased recognition and application of 
antiferroelectrics as functional materials. 

An antiferroelectric [TJ [8] is like a ferroelectric in that 
its structure is obtained through distortion of a nonpolar 
high-symmetry reference phase; for ferroelectrics, the dis- 
tortion is polar, while for antiferroelectrics it is nonpolar. 
However, not all nonpolar phases thus obtained are an- 
tiferroelectric: in addition, there must be an alternative 
low-energy ferroelectric phase obtained by a polar distor- 
tion of the same high-symmetry reference structure, and 
an applied electric field must induce a first-order transi- 
tion from the antiferroelectric phase to this ferroelectric 
phase, producing a characteristic P-E double-hysteresis 
loop |9j [10] . This behavior in applied electric field is the 
origin of the functional properties of antiferroelectrics. 
For systems in which there is a difference in lattice pa- 
rameters between the AFE and FE phases, the electric- 
field-induced transition produces a large nonlinear strain 
response, useful for transducer applications [11 . The en- 



tropy change between the two phases similarly can pro- 
duce a high effective electrocaloric response, useful for 
solid-state cooling [12]. Finally, the shape of the P-E 
hysteresis loop allows storage of electrical energy of suf- 
ficient magnitude to attract interest for energy storage 
applications p~3j [14] . 

Most of the research on antiferroelectric oxides has 
focused on two classes of materials, both with rather 
complex structures. The Pbam structure of the proto- 
typical perovskite antiferroelectric PbZrOs [9] is obtained 
through a cell-quadrupling antipolar Pb displacement 
mode combined with four additional modes that further 
double the unit cell [15]; the Pbcm structure of the an- 
tiferroelectric phases of NaNbOs and AgNbOa are simi- 
larly complex [16]. Chemical substitution into the end- 
point compounds allows tuning of the critical temper- 
ature, the critical field, the electric-field induced strain 
and polarization, and other functional properties [THE]- 
Chemical substitution is also seen to induce transitions to 
a distinct tetragonal antiferroelectric phase, for example 
in Pbi_ ;E Sr ;E (Zri_ 2y Ti 2/ )03 Jl , and to related ferroelec- 
tric phases, for example in PbZri-^Ti^Oa, also obtained 
as distortions of the ideal perovskite structure. 

Evidence is beginning to emerge that antiferroelectric- 
ity in inorganic materials extends far beyond the systems 
which have been the focus of the literature to date. In 
particular, recent observations of orthorhombic Pnma 
Sm-doped BiFe0 3 [E] and thin-film BiCr0 3 [19 show 
the double hysteresis loops characteristic of antiferro- 
electricity. This suggests that many additional nonpolar 
phases should be considered as candidates for previously- 
unrecognized antiferroelectricity. Our strategy to search 
for new antiferroelectric materials is to identify a class of 
materials with a structure type that is obtained from an 
antipolar distortion of a high-symmetry reference struc- 
ture which is also related through a polar distortion to 
a structure type of known ferroelectrics. By screen- 
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ing a large number of both reported and hypothetical 
compounds, we can find compositions that are insulat- 
ing and locally stable both in the antiferroelectric struc- 
ture and in the related ferroelectric structure; of these, 
the compounds for which the antiferroelectric structure 
is slightly lower in energy than the ferroelectric struc- 
ture are promising candidates for new antiferroelectrics 
for targeted experimental investigation. Through this 
search, we expect to find new AFE structure types which 
are significantly simpler than the Pbam structure, even 
perhaps including realizations of the two-sublattice Kit- 
tel model [10] and q=0 antiferroelectrics [20]. For tech- 
nological applications, new materials offer the possibil- 
ity of better performance via reduced hysteresis, larger 
changes at the electric-field-induced transition, operation 
at higher or lower temperatures, better integration with 
other materials based on structural or chemical compat- 
ibility, elimination of toxic elements, and/or a band gap 
suitable for photoactive applications. 

ABC compounds in the Pnma MgSrSi structure 
type [21 (previously referred to as the anti-PbCl2 struc- 
ture type [22] ) are a promising target class. The MgSrSi 
structure, shown in Fig. [I] is obtained by a nonpolar 
distortion of the high-symmetry P63 / mmc ZrBeSi struc- 
ture. The distortion can be decomposed into two modes: 
an M2 mode (antipolar displacements of ABC along c, 
and antipolar displacements of A along b) that breaks the 
Pfts/mmc symmetry to Pnma, followed by a Tl~ mode 
(shifting of BC layers along b) that does not break any 
additional symmetry. A polar distortion specified by a 
r^" mode relates the high-symmetry P^jmmc structure 
to the LiGaGe-type Purine ferroelectric compounds iden- 
tified previously [5] ; in addition to being insulating, these 
compounds can have spontaneous polarization compara- 
ble to that of BaTiOs. However, little has been reported 
about the band structure of MgSrSi-type compounds or 
about their response to applied electric fields [23] [24] . 

In this paper, we use first-principles methods to es- 
tablish a new class of antiferroelectrics in the MgSrSi 
structure type and to identify promising candidate ma- 
terials for further investigation. Specifically, we compute 
the structural parameters, nonpolar distortions, band 
gaps, and AFE-FE energy differences for a search set 
comprising 37 reported and 33 as-yet-hypothetical ABC 
compounds in the MgSrSi structure type. We iden- 
tify 11 combinations for which both the FE and AFE 
phases are insulating and have an energy difference be- 
low 0.2 eV. For all insulating combinations studied, we 
find that the band gaps are in the semiconducting range; 
the lower band gaps could be useful for photoactive appli- 
cations [25 27 . These candidate antiferroelectrics offer 
promise for experimental investigation and for the future 
development of new high-performance materials for prac- 
tical applications. 

First principles computations were performed with the 
ABINIT package [28] . The local density approximation 
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FIG. 1: Left: planar buckling distortions of the polar PQ^mc 
LiGaGe structure relative to the P63/mmc high symmetry 
ZrBeSi structure type (center), as described in Ref. [5 . Right: 
antipolar distortions of the Pnma MgSrSi structure type rel- 
ative to the high symmetry ZrBeSi structure type (center), 
showing antipolar buckling of the planes formed by atoms 
BC at Wyckoff positions 2b (dark blue) and 2b' (gold) and 
antipolar displacements of the stuffing atoms A at 2a (green) . 

(LDA) and a 6x6x6 Monkhorst-Pack sampling of the 
Brillouin zone [29] were used for all calculations, ex- 
cept for the Berry phase polarization [30j [31] calcula- 
tions, for which an 8x8x8 grid was used. All atoms 
were represented by norm-conserving, optimized [32] . de- 
signed nonlocal [33] pseudopotentials, generated with the 
OPIUM code [34] . All calculations were performed with a 
plane wave cutoff of 50 Ry. In addition, the QUANTUM 
ESPRESSO [35] package was used to perform nudged 
elastic band calculations [S3 [37]. WANNIER90 [38] was 
used to generate maximally localized Wannier functions 
(MLWF) [39]. 

Our search set of candidate MgSrSi-type antiferro- 
electrics consists of ABC combinations drawn from the 
Inorganic Crystal Structural Database (ICSD) [40 and 
two recent high-throughput searches. The ICSD in- 
cludes 37 previously synthesized ABC compounds in the 
MgSrSi structure type that do not contain an /-block 
element and which have a total of 8 or 18 valence elec- 
trons, which promotes band gap formation. In addition, 
we included 13 as- yet unsynthesized compounds from the 
high-throughput study of Zhang et al. [3 that are pre- 
dicted to have the desired Pnma structure. Finally, we 
included 20 compounds from the high-throughput study 
of Bennett et al. [5 J that are predicted to be insulating 
and at least locally stable in the polar LiGaGe structure 
type, bringing our total search set to 70 compounds. We 
can classify these combinations into the following groups: 

I- I-VI (NaKSe), I-II-V (LiCaBi), I-XII-V (NaCdAs), XI- 

II- V (CuMgP), II-II-IV (SrCaGe), III-X-V (ScNiP), III- 
XI-IV (ScCuSi), and IV-X-IV (ZrPtSi). 

For each ABC combination in our search set, we op- 
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timize the structural parameters for each of the three 
structural variants ABC, ABC, and ABC, where the 
underscore indicates the stuffing atom (Fig. [T]), for each 
of our three structure types: MgSrSi (Prima), LiGaGe 
(PQ^mc), and ZrBeSi (P^/mmc). The computed struc- 
tural parameters (Supplemental Table 1) generally show 
good agreement with experimental values, with the un- 
derestimate of lattice constants characteristic of LDA cal- 
culations, about 1-3 % for a and as large as 3-4 % for 
c[44j. While most of our calculations agree with previ- 
ous experimental and theoretical determinations of the 
ground state structure, there are a few minor discrepan- 
cies. We find that NaCdSb prefers the P6smc structure 
by 34 meV/f.u. over the experimentally reported Prima 
structure. In addition, we find that KSrBi and RbBaAs 
prefer the high-symmetry P6s/mmc structure, in con- 
trast to the results of Zhang et al.[3], which we attribute 
to minor methodological differences [45] . 

Our most promising candidate antiferroelectrics, 
shown in Table [l| fulfill the following criteria: a) the 
ground state is the antipolar Prima MgSrSi structure 
type, b) the polar P6smc LiGaGe structure type is 
metastable, and c) both phases are insulating. In Table [i] 
we present the energy differences between the Prima and 
P6smc structures (AE), the energy difference between 
the P§zmc and P6s/mmc (AEsw), the polarization of 
the LiGaGe structure, the mode decomposition of the an- 
tipolar distortions, and difference in volume between the 
P6smc and Prima structures. Of this set, seven com- 
binations are known materials, highlighted in boldface 
in Table [TJ The calculated band gaps of these materi- 
als cover the entire semiconducting range, from 0.04 eV 
(BaCaSi) to 2.91 eV (NaLiTe), with changes of 0.01 to 
0.59 eV upon switching (see Supplemental Tables). 

We find eleven insulating compounds with energy dif- 
ferences between the polar and antipolar phases below 
0.2 eV/f.u. (see Table 0. Of this set, six combinations, 
namely LiCaBi, LiCaSb, KNaTe, NaLiTe, KNaS, and 
KNaSe, are known materials. Despite the small energy 
differences, many of our candidate materials still require 
a large critical field to stabilize the polar phase. How- 
ever, the structural versatility of ABC intermetallic com- 
pounds should make it experimentally feasible to reduce 
the critical field of many of these materials via formation 
of solid solutions of chemically similar compounds with 
different structures. For instance, we find that LiBeP has 
an antipolar Prima ground state, while LiBeAs has a po- 
lar P6smc ground state. A solid solution of these two 
materials should result in a material in which it is possi- 
ble to chemically tune towards the first-order transition 
between preferred structures, thus reducing the critical 
field. 

In addition to our most promising candidates, in Ta- 
ble |TT] we present data for seven compounds that have 
a Prima ground state as well as a stable polar distor- 
tion, but for which only one of the two phases is in- 



ABC 


AE S w 


P 


r 5 + 


M- 


AE 


AV 


AV/V 




(meV) 


(C/m 2 ) 






(meV) 


(A 3 ) 


(%) 


UBeP 


119 


0.85 


0.26 


1.26 


18 


1.23 


3.7 


MgLiP 


20 


0.38 


0.41 


1.28 


230 


0.49 


1.1 


MgLiAs 


30 


0.39 


0.42 


1.35 


207 


0.12 


0.2 


CaLiSb 


7 


0.18 


0.19 


1.17 


79 


0.09 


0.1 


CaLiBi 


7 


0.19 


0.20 


1.16 


80 


-0.10 


-0.1 


NaMgP 


102 


0.49 


0.31 


1.43 


275 


1.88 


3.2 


NaMgAs 


114 


0.48 


0.34 


1.44 


232 


1.65 


2.5 


NaMgSb 


146 


0.43 


0.40 


1.51 


154 


0.93 


1.2 


NaMgBi 


127 


0.42 


0.44 


1.51 


143 


1.14 


1.4 


KMgSb 


41 


0.40 


0.33 


1.77 


254 


2.17 


2.6 


KMgBi 


73 


0.31 


0.34 


1.78 


227 


1.65 


1.8 


NaZnSb 


81 


0.49 


0.18 


1.48 


42 


2.68 


4.1 


NaLiTe 


23 


0.20 


0.30 


1.41 


114 


0.96 


1.4 


KNaS 


12 


0.17 


0.32 


1.36 


149 


1.60 


2.3 


KNaSe 


13 


0.15 


0.33 


1.38 


131 


2.09 


2.7 


KNaTe 


14 


0.13 


0.38 


1.42 


96 


2.81 


2.9 


NaAgSe 


51 


0.65 


0.62 


1.63 


134 


0.11 


0.2 


BaCaSi 


11 


0.34 


0.35 


1.62 


310 


0.37 


0.4 



TABLE I: Energy difference between Purine and PQ^/mmc 
phases (A^sw), polarization P, amplitude of the microscopic 
distortions from ZrBeSi to MgSrSi structure types (r^ and 
M 2 ~), energy difference between AFE and FE states (AE), 
and change in volume (AV), also given as a percentage, for 
the eighteen compounds in the search set that are insulating in 
both the antipolar and polar states. Combinations reported 
as Prima in the ICSD are shown in boldface. Energies are 
reported in meV per formula unit. Band gaps are reported in 
the supplemental material. 
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TABLE II: Same as Table 1, but for compounds in the search 
set for which either the Prima state or the P6smc state is 
metallic in our DFT calculation. A dash in the polarization 
(P) column indicates that the PQ^mc state is metallic. 



sulating. This caveat would normally disqualify these 
compounds as candidate antiferroelectrics, but it is well 
known that DFT-LDA tends to significantly underesti- 
mate band gaps, and can even predict small-gap semicon- 
ductors (e.g., Ge) to be metallic [41]. It is possible that 
these seven compounds could display antiferroelectric be- 
havior, and therefore merit experimental investigation, 
especially since four are known materials (highlighted in 
boldface in Table [ll|. 
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FIG. 2: Minimum-energy path for switching between the non- 
polar Pnma and polar PQ^mc phases of NaCdAs, as calcu- 
lated by the NEB method, (a-c): isosurfaces of an inter- 
planar bonding state (MLWF), which reverses orientation 
during the switching process. Lower panel: graph of energy 
vs reaction coordinate, showing that our proposed reaction 
pathway from Pnma to PQ^mc has a barrier of 0.11 eV. 
This is in comparison to a path that would have the high- 
symmetry PQz/mmc structure as our intermediate, which 
would be 0.19 eV. 



In order to investigate the switching path between the 
antipolar and polar states, we use nudged elastic band 
calculations (with the volume fixed to the Pnma cell) to 
identify the lowest-energy path between the two phases, 
as shown in Fig. [2] While we expect that the actual 
experimental switching path will depend on complex ki- 
netic processes like domain nucleation and domain-wall 
motion, our calculations of homogeneous switching paths 
provide an estimate of the energy scales involved. For 
our example of NaCdAs in Fig. [2j we find that the low- 
energy switching path from Pnma to P6smc corresponds 
primarily to half of the Cd atoms moving from below the 
As layer to above it, while the Na atoms rearrange to 
avoid the Cd. For NaCdAs, this first-order transition 
path reduces the barrier to switching by more than 50% 
as compared to switching through the high-symmetry 
P6s/mmc structure, and makes it more likely that the 
material will return to the ground state when the exter- 
nal field is switched off. 

We use maximally localized Wannier functions to an- 
alyze differences in bonding between the P^mc and 
Pnma structures and to study how the bonding changes 
when the material is switched. While we refer to the 
previously synthesized NaCdAs as our example, we find 
that our description of bonding is similar in other combi- 
nations. The bonding orbitals of the P6s mc structure 
are sp 3 -like orbitals centered near As, the most elec- 
tronegative element, with all bonds oriented towards the 



four neighboring Cd atoms (see Fig. [2|. The bonds of 
the Pnma structure are also sp 3 -like and centered on 
As; however, a transition between the two structures 
requires that half of the sp 3 bonds reorient. This pro- 
ceeds by breaking an inter-layer Cd-As bond, with this 
orbital forming a p z -\ike non-bonding state in the transi- 
tion structure, and then reforming the bond pointing in 
the other direction (see the orbital in Fig. [2|. In addi- 
tion, the in-plane sp 3 bonds reorient by passing through 
an sp 2 -\ike intermediate state. The similar bonding in 
the P6smc and Pnma structures is consistent with their 
competitive energies. However, we note that in the 
Pnma structure, only half of the sp 3 bonds are oriented 
directly between Cd and As, which accounts for the sig- 
nificant changes in band gap between the two structure 
types. 

Our proposed antiferroelectrics have many properties 
that could make them better for applications than ex- 
isting materials. For instance, most of our candidates 
have a change in volume of 1-4% upon switching (see 
Table [i]), which is much larger than the change observed 
in most piezoelectrics (0.1-0.2%) and current antiferro- 
electrics (< 0.9%) [TT] , Very large volume changes could 
make these materials ideal for a variety of transducer 
applications in which harnessing a large non-linear re- 
sponse is necessary. Another application would be to use 
these antiferroelectrics as a substrate which could apply 
a reversible strain to thin films grown on top. We also 
note that in contrast to perovskite oxides, the stability 
of the polar and antipolar distortions in these materials 
is relatively insensitive to strain [5], these materials have 
no competing distortion patterns, and they are strongly 
anisotropic, a combination of favorable properties which 
could make antiferroelectric ABC useful in applications 
that require a robust and reversible response to exter- 
nal fields, such as high-energy storage capacitance and 
electrocaloric refrigeration. 

In conclusion, we have used first-principles methods to 
establish a new class of antiferroelectrics in the MgSrSi 
structure type and to identify promising candidate ma- 
terials for further investigation. Through targeted syn- 
thesis, MgSrSi-type compounds could potentially be de- 
veloped as a valuable class of functional antiferroelectric 
materials. This is a specific application of a larger-scale 
strategy to identify new functional materials by targeting 
insulating compounds not previously recognized as func- 
tional materials and tuning the composition and other 
control parameters, such as epitaxial strain, and/or mod- 
ifying the structure by intercalation of atoms. The iden- 
tification of antiferroelectricity in classes of materials in 
which it was previously unrecognized offers the possibil- 
ity of optimizing properties and combining polarization 
with other functional properties, including magnetism, 
to produce multifunctional behavior of fundamental sci- 
entific interest and for groundbreaking technological ap- 
plications. 
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